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ABSTRACT 

Context. During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important 
feedback processes. Both processes reflect the accretion onto the forming star and affect subsequent collapse or fragmentation of the 
cloud. 

Aims. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources in a consistent manner. 
The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations 
and evolutionary trends. 

Methods. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the 
CHAMP^ instrument on the Atacama Pathfinder Experiment (^^CO and 6-5; Eup^lOO K), and the HARP-B instrument on the 
James Clerk Maxwell Telescope (^^CO and 3-2; Eup~30 K). The maps have high spatial resolution, particularly the CO 6-5 
maps taken with a 9" beam, resolving the morphology of the outflows. The maps are used to determine outflow parameters and the 
results are compared with higher-/ CO lines obtained with Herschel. Envelope models are used to quantify the amount of UV-heated 
gas and its temperature from ^^CO 6-5 observations. 

Results. All sources in our sample show outflow activity, with the spatial extent decreasing from the Class 0 to the Class I stage. 
Consistent with previous studies, the outflow force, Fco, is larger for Class 0 sources than for Class I sources, even if their luminosities 
are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences 
the surrounding cloud material directly. Comparison of the CO 6-5 results with HIFI H 2 O and PACS high-/ CO lines, both tracing 
currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link 
does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas 
CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing 
the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. 
The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as 
important on small scales. 

Key words. Astrochemistry — Stars: formation — Stars: protostars — ISM: molecules — Techniques: spectroscopic 


1. Introduction 


During the early phases of star-formation, material surrounding 
the newly forming star accretes onto the protostar. At the same 
time, winds or jets are launched at supersonic speeds from the 
star-disk system, which sweep up surrounding envelope material 
in large bipolar outflows. The material is accelerated and pushed 
to distances of several tens of thousands of AU, and these out¬ 
flows play a pivotal role in the physics and chemistry of the star¬ 
forming cores ( [Snell et aL]|1980[ [Goldsmith et al.|p^84} [Lada 


1987l[Greene etal.|1994l 

[Bachiller & Tafalla| 1999 [[Arce & Sar-[ 

gent[f2006[ Tafalla et al.| 

[2013[). The youngest protostars have 


highly collimated outflows driven by jets, whereas at later stages 
wide-angle winds drive less collimated outflows. However, there 
is still not a general consensus to explain the launching mecha¬ 
nisms and nature of these outflows ([Arce et al.|2007l[Erank et al.[ 

[MT^ . 

The goal of this paper is to investigate how the outflow ac¬ 
tivity varies with evolution and how this compares with other 
measures of the accretion processes for low-mass sources. The 
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outflows reflect the integrated activity over the entire lifetime of 
the protostar, which could be the result of multiple accretion and 
ejection events. It is important to distinguish this probe from the 
current accretion rate, as reflected for example in the luminosity 
of the source, in order to understand the accretion history. The 
well-known “luminosity problem” in low-mass star-formation 
indicates tha t protostars are underluminous compared to theoret¬ 
ical models ([Kenyon et al.|1990t Evans et al.|2Q09t [Enoch et al 


2009[ [Dunham et aT|20l6[ |2013| ). One of the possible resolu¬ 

tions to this problem is that of “episodic accretion”, in which the 
star builds up through short bursts of rapid accretion over long 
periods of time rather than continuous steady-state accretion. An 
accurate and consistent quantiflcation of outflow properties, such 
as the outflow force and mass, is essential for addressing this 
problem. 

Outflows have been observed in CO emission in the last 
few decades towards many sources, but those observations were 
mainly done via lower-/ CO rotational transitions (/u<3), which 
probe colder swept-up or entrained gas (r~50-10Q K) (e.g.. 


Bachiller et al.p^^ [Blake et al.||1995| [Bontemps et al.||1996[ 

Tafalla et al.|2000[ [Curtis et al.|2010[ and many others). One of 

the most important parameters that is used for the evolutionary 
studies of star formation is the “outflow force”, which is known 
as the strength of an outflow and deflned similar to any r“^-type 
force. These studies conclude that the outflow force correlates 
well with bolometric luminosity, Lboi, a correlation which holds 
over several orders of magnitude. Furthermore, the outflow force 
from Class 0 sources is stronger than for Class I sources, indicat¬ 
ing an evolutionary trend. The correlations, however, often show 
some degree of scatter, typically more than an order of magni¬ 
tude in Fco for any value of Lboi- Some of the uncertainties in 
these studies include the opacity in the line wings, the adopted 
inclination angle and cloud contamination at low outflow veloci¬ 
ties (e.g., [van der Marel et al.[201^ . Comparison with other out¬ 
flow tracers such as water recently observed with the Herschel 
Space Observatory is further complicated because the various 
studies use different analysis methods to derive outflow param¬ 
eters from low-/ CO maps. One of the goals of this paper is to 
provide a consistent set of outflow parameters determined by the 
same method using data from the same telescopes for compari¬ 
son with the Herschel lines. 

Recently, the importance of radiative feedback from low- 
mass protostars on all scales of star formation has been acknowl¬ 
edged. On cloud scales (>10"^ AU) the feedback sets the effi¬ 
ciency at which cores fragment from the cloud and form stars 
( [Qffner et al.|2009 [2010[ [Hansen et al.|2012[ ) because the Jeans 
length scales as . Simulations including radiative feedback 
and radiative transfer reproduce the observed initial mass func¬ 
tion (IMF) better than models without these effects included 
( [Qffner et al.[[2009[ ). On the scales of individual cores (<3000 
AU), the radiative feedback suppresses the fragmentation into 
multiple systems a nd serves to stabilize the protostellar disk 
( [Qffner et al.|2010[ ). Thus, quantifying observationally the tem¬ 
perature changes as a function of position from the protostar 
are important steps toward more accurate models of star for¬ 
mation. The first observational evidence of heating of the gas 
around low-mass protostars on scales of -1000 AU by UV radi¬ 
ation escaping through the outflow cavities dates back to [Spaans[ 
et al. (1995) based on strong narrow ^^CO 6-5 lines, and has 
since been demonstrated and quantified for a few more sources 
by[van Kempen et^ ( [20091^ ; [Yildiz et al.[ ( [201^ ; [Visser et al.[ 
( [2Q12[ ). We note that this UV-heated gas is warm gas with tem¬ 
peratures higher than that of the dust, and is thus in excess of 
warm material in the envelope that has been heated by the proto- 

Article number, page 2 of[^ 


stellar luminosity, where the gas temperature is equal to the dust 
temperature. Although UV heating toward photo-dissociation re¬ 
gions (PDRs) is readily traced by emission from polycyclic aro¬ 
matic hydrocarbons (PAHs), the PAH abundance toward embed¬ 
ded protostars is too low for them to be used as a tool in this con¬ 
text ( [Geers et ar][2009[ ). Here we investigate the importance of 
radiative feedback for a much larger sample of low-mass sources 
and compare the gas temperatures and involved mass with that 
of the outflows. 

Tracing warm gas (r>30 K) in the envelope or in the sur¬ 
roundings requires observations of higher-/ transitions of CO, 
e.g., /u>5, for which ground-based telescopes demand excellent 
weather conditions on dry observing sites. The CHAMP^ instru¬ 
ment, mounted on the Atacama Pathfinder Experiment (APEX) 
telescope is ideally suited to observe higher-/ CO transitions 
and efficiently map extended sources. The broad line wings 
of CO 6-5 (£’u/k=115K) suffer less from opacity effects than 
CO 3-2 (£’u/k=33K) ( [van Kempen et al. [2009a[ [Yildiz et al. 
[2012[ ). Moreover, the ambient cloud contribution is smaller for 
these higher-/ transitions, except close to the source position, 
where the dense protostellar envelope may still contribute. Even 
higher-/ CO lines up to /u~50 were routinely observed with 


the shocked gas in the Herschel beam (I 

Jerczeg et al. 2012[ 

Goicoechea et al.[[2012[ 

Benedettini et a . 

[2012[ [Manoj et al. 

2013t|Greenetal.|2013l| 

Nisini et al.|2013[ 

Karska et al.|2013[). 


This currently shocked gas is different from that observed in 
low-/ CO transitions, as is evident from their different spatial 
distributions ( [Tafalla et al.|2013[|Santangelo et al.|20T3] ). 

In this paper, we present an APEX-CHAMP^ survey of 26 
low-mass young stellar objects (YSOs), which were mapped in 
CO /=6-5 and isotopologues in order to trace their outflow ac- 
tivity, following f^n Kempen et al.[ p009a|b[ ) and [Yildiz et al.[ 


( [2012[ ), papers I, II and III in this series, on individual or more 
limited samples of sources. These data complement our earlier 
surveys at lower frequency of CO and other molecules with the 


James Clerk Maxwell Telescope (JCMT) and APEX (e.g., [j0r- 
[gensen et al.[[2002[[2QQ4[[\^ Kempen et al.[[2009c[ ). The same 
sources are covered in the Herschel key pr oject, “Water in star- 
forming regions with HerscheF (WISH; [van Dishoeck et aT 
[201 1[ ), which has observed H 2 O and selected high-/ CO lines 
with HIFI and PACS instruments. Many of the sources are also 
included in the “Dust, Ice and Gas in Time” program (DIGIT; PI: 
N. Evans; [Green et al.[2013[ ), which has obtained full PACS spec¬ 
tral scans. The results obtained from the ^^CO maps are comple¬ 
mented by ^^CO 6-5 data of the same sources, with the narrower 
^^CO 6-5 lines probing the UV photon-heated gas. 

The YSOs in our sample cover both the deeply em bedded 
Class 0 stage as well as the less embedded Class I stage ( [Andre 
[et al.pOOOt [Robitaille et al.[[2006[ ). Physical models of the dust 
temperature and density stru cture of the envelopes h ave been 
developed for all sources by [Kristensen et al.[ ( [2012[ ) through 
spherically symmetric radiative transfer models of the contin¬ 
uum emission. The full data set covering many sources, together 
with the envelope models, allows us to address important char¬ 
acteristics of YSOs through the evolution from Class 0 to Class I 
in a more consistent manner. These characteristics can be in¬ 
ferred from their different morphologies, outflow forces, enve¬ 
lope masses, etc. and eventually be compared with evolutionary 
m odels. The study pre sented here is also complementary to that 
of [Yildiz et al.[ ( [2013] ), where only the source position was stud¬ 
ied with spectrally resolved CO line profiles from /=2-l to 10-9 
(£’up~300 K), and trends with evolution were examined. 
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Fig. 1. Envelope mass,Menv, vs bolometric luminosity, Lboi, for the sur¬ 
veyed sources. Red diamonds and blue squares indicate Class 0 and 
Class I sources, respectively. 


The outline of the paper is as follows. In Section the ob¬ 
servations and the telescopes where the data have been obtained 
are described. In Section ^ physical parameters obtained from 
molecular outflows are given and the UV heated gas component 
is identifled. In Sect. these results are discussed, and conclu¬ 
sions from this work are presented in Sect. 


2. Sample and observations 


2.1. Sample 

The sample selection criteria with the coordinates and other ba¬ 


sic information of the source list are presente d in [van Dishoeck 
et al.| ( |MTT] ) with updates in |Kristensen et al.|(|20i2 ), and is the 


same as the sample presented in |Yildiz et al. ( |2Q13| ). It consists 
of 15 Class 0 and 11 Class I embedded protostellar sources lo¬ 
cated in the Perseus, Ophiuchus, Taurus, Chamaeleon, and Ser¬ 
pens molecular clouds. The average distance is 200 pc, with a 
maximum distance of 450 pc. 

Figure presents the envelope mass (Menv) as a function of 
bolometric luminosity (Lboi) for all sources. The parameters are 
taken from the continuum radiative transfer modeling by |Kris- 


tensen etaL] ( |2012| ) based on fits of the spectral energy distri¬ 
butions (SEDs) including new Herschel-FACS fluxes, as well 
as the spatial extent of the envelopes observed at submillime¬ 
ter wavelengths. The envelope mass is measured either at the 
T'dust^lO K radius or at the ^=10"^ cm“^ radius, depending on 
which is smaller. Class 0 and Class I sources are well separated 
in the diagram, with the Class 0 sources having higher envelope 
masses. This type of correlation diagram has been put forward 
by |Saraceno et al.| ( [l996| ) and subsequently used as an evolution¬ 
ary diagram for embedded YSOs with lower envelope masses 
representing lat er stages (e.g.,[Bontemps et al.|1996} [Hogerhei^ 


[jde et al.||1998 [Hatchell et al.||2007| ). In our sample, envelope 

masses range from 0.04 M© (Elias 29) to 16 M© (Ser-SMMl) 
and the luminosities range from 0.8 L© (CedllOIRS4) to 35.7 
L© (NGC 1333-IRAS 2A). The large range of masses and lumi¬ 
nosities makes the sample well suited for studying trends with 
various source parameters. The range of luminosities studied is 
similar to that of |Bontemps et al.| ( [1996| ), -0.5 to 15 L©, but our 


sample is more weighted toward higher luminosities and earlier 
stages. 


2.2. Observations 


Molecular line observations of CO in the /=6-5 transitions were 
done with the 12-m submillimeter A tacama Pathflnder Experi¬ 
ment, APE^ ( [Gusten et al. 2008 ) at Llano de Chajnantor in 
Chile, whereas the J=3-2 transition was primarily observed at 
the 15-m James Clerk Maxwell Telescope, JCMT0 at Mauna 
Kea, Hawaii. 

APEX: ^^CO and ^^CO 6-5 maps of the survey were ob¬ 
tained with the CHAMP^ instrument on APEX between June 
2007 and September 2012. The CHAMP^ instrument consists of 
two heterodyne receiver arrays, each with seven pixel detector 
elements for simultaneous operations in the 620-720 GHz and 
780-950 G Hz frequency ranges ( [Kasemann et al.|2006[ |Gusten 
et al.||2008 1. The observational procedures are explained in de 


tail in [van Kempen et aL] ( |2009a|b|c| ) and [Yildiz et aL] ( |2012 ). 
Simultaneous observations were done with the following set¬ 
tings of the lower and higher frequency bands: ^^CO 6-5 with 
^^CO 7-6; ^^CO 6-5 with [Ci] 2-1. ^^CO maps cover the entire 
outflow extent with a few exceptions (L1527, CedllOIRS4, and 
L1551-IRS5), whereas ^^CO maps cover only a ~100"xl00" 
region around the central source position. LI 157 is part of the 
WISH survey, but because it is not accessible from APEX (dec 
= +68°), no CO 6-5 data are presented. 

The APEX beam size is -9" (-1800 AU for a source 
at 200 pc) at 691 GHz. The observations were done using 
position-switching toward an emission-free reference position. 
The CHAMP^ instrument uses the Past Pourier Transform Spec¬ 
trometer (PETS) backend ( [Klein et al.|2006| ) for all seven pixels 
with a resolution of 0.183 MHz (0.079 km s“^ at 691 GHz). The 
rms at the source position is listed in [Yildiz et al.| ( |20l'3] ) for the 
CO 6-5 and ^^CO 6-5 observations and is typically 0.3-0.5 K 
for the former and 0.1-0.3 K for the latter, both in 0.2 km s“^ 
channels. The rms increases near the map edges where the effec¬ 
tive integration time per beam was significantly smaller than in 
the central parts; near the edges the rms may be twice as high. 
Apart from the high-/ CO observations, some of the 3-2 line 
observations were also conducted with APEX for a few southern 
sources, e.g., DK Cha, CedllOIRS4, and HH46 
[et al.|2009c[ ). 

JCMT: Pully sampled jiggle maps of ^^CO and ^^CO 3-2 
were obtained using the HARP-B instrument mounted on the 
JCMT. HARP-B consists of 16 SIS detectors with 4x4 pixel el¬ 
ements of 15" each at 30" separation. Most of the maps were 
obtained through our own dedicated proposals, with a subset ob¬ 
tained from the JCMT public archiv^ 

The data were acquired on the antenna temperature 
scale and were converted to main-beam brightness temperatures 
Tmb=T^/^mb using the beam efficiencies (i/mb)- The CHAMP^ 


van Kempen 


^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between the 
Max-Planck-Institut fur Radioastronomie, the European Southern Ob¬ 
servatory, and the Onsala Space Observatory. 

^ The JCMT is operated by The Joint Astronomy Centre on behalf of 
the Science and Technology Eacilities Council of the United Kingdom, 
the National Research Council of Canada, and (until 31 March 2013) 
the Netherlands Organisation for Scientific Research. 

^ This research used the facilities of the Canadian Astronomy Data 
Centre operated by the National Research Council of Canada with the 
support of the Canadian Space Agency. 
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Fig. 2. Overview of the outflows traced by the 6-5 observations with the APEX-CHAMP^ instrument. Contour levels are given in Table A.l 


and the source is located at (0,0) in each map, with the exception of the maps of NGC 1333-IRAS 4A and IRAS 4B, and Ser-SMM3 and Ser- 
SMM4, which are located in the same maps and centered on NGC 1333-IRAS 4A and Ser-SMM3, respectively. The circle in each plot c orresp onds 
to a region of 5000 AU radius at the distance of each source. Velocity ranges over which the integration was done are provided in Table] A.l] 


beam efficiencies were taken from the CHAMP^ websit^ and 
forward efficiencies are 0.95 in all observations. The various 
beam efficiencies are all given in |Yildiz ctH] ( |2Q13[ their Ap¬ 
pendix C) and are typically -0.5. The JCMT beam efficiencies 
were taken from the JCMT Efficiencies Databas^ and 0.63 is 
used for all HARP-B observations. Calibration errors are esti¬ 
mated to be -20% for both telescopes. Typical rms noise levels 
of the J=3-2 data are from 0.05 K to 0.1 K in 0.2 km s“^ chan¬ 
nels. 

For the data reduction and analysis, the “Continuum and 
Line Analysis Single Dish Software”, CLASS program, which is 
part of the GILDAS softwar^ is used. In particular, linear base¬ 
lines were subtracted from all spectra. ^^CO and ^^CO 6-5 and 
3-2 line profiles of the central source positions of all the sources 
in the sample are presented in |Yildiz et al.| ( |201? l. 


2.3. CO maps 

All spectra are binned to a 0.5 km s“^ velocity resolution for an¬ 
alyzing the outflows. The intensities of the blue and red outflow 
lobes are calculated by integrating the blue and red emission in 
each of the spectra separately, where the integration limits are 

http://www3.mpifr-bonn.mpg.de/div/submintech/ 
heterodyne/champplus/champ_efficiencies.29-ll-13.html 
^ http: //WWW. jach.hawaii.edu/JCMT/spectral_line/ 
Standards/eff_web.html 
” http: //WWW .iram.fr/IRAMFR/GILDAS 
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carefully selected for each source by using the 0.2 km s~^ reso- 
lutionCO 3-2 or 6-5 spectra if the former is not available (see 
Fig. |A.l| ). First, the inner velocity limit, Vin, closest to the source 
velocity is determined by selecting a spatial region not associ¬ 
ated with the outflow. The ^^COspectra in this region are av¬ 
eraged to determine the narrow line emission coming from the 
envelope and surrounding cloud, and Vin is estimated from the 
width of the quiescent emission (see Fig. |A.l| in the Appendix). 
Second, the outer velocity limits Vout are determined from the 
highest S /N spectrum inside each of the blue and red outflow 
lobes. The outer velocity limits are selected as the velocity where 
the emission in the spectrum goes down to the Icr limit for the 
first time. It therefore excludes extremely high velocity or ‘bul¬ 
let’ emission which is seen for a few sources. The blue- and 
red-shifted integrated intensity is measured by integrating over 
these velocity limits across the entire map, but excluding any 
extremely high velocity (EHV) or “bullet” emission. 


2.3.1. Outflow velocity 

The maximum outflow velocity, Ymax is defined as IKut-kLSRL 
the total velocity extent measured relative to the source velocity. 
In order to estimate Vmax, representative spectra from the blue 
and red outflow lobes observed in CO 3-2 are selected sepa¬ 
rately, and Vout is measured as described above. The differences 
between the velocity, where the emission reaches Icr level (Vout) 
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Fig. 3. Overview of the entire set of outflows traced by the 3-2 observations with the JCMT and APEX. Contour levels are given in Table | A. l] 
and the source is located at (0,0) in each map, with the exception of the maps of NGC 1333-IRAS 4A and NGC 1333-IRAS 4B, and Ser-SMM3 
and Ser-SMM4, which are located in the same maps and centered on NGC 1333-IRAS 4A and Ser-SMM3, respectively. The circle in each plot 
corre spond s to a region of 5000 AU radius at the distance of each source. Velocity ranges over which the integration was done are provided in 
TableEI] 


with Vlsr are taken as the global V^ax values for the correspond¬ 
ing blue and red-shifted lobes ( [Cabrit & Bertout|1992| ). 

Two issues arise when determi ning V^ax (e.g., |van der Marel] 
|et al.||2013l [Dunham et al.||2Ql^ : First, V^ax is a function of 
the rms noise level and generally decreases with increasing rms. 
For noisy data, Vmax may be underestimated compared to its true 
value. For this reason, the 3-2 lines are chosen to determine V^ax 
because of their higher S/N than the 6-5 lines. Second, if the 
outflow lobes are inclined, Fmax suffers from projection effects. 
Both effects will increase the value of Vniax if properly taken into 
account. 

Concerning the second issue, the inclination is difficult to es¬ 
timate from these data alone; proper-motion studies along with 
radial velocities are required to obtain an accurate estimate of 
the inclination. Alternatively, the velocity structure may be mod¬ 
eled assuming some distribution of material, e.g., a wind-driven 
shell with a Hubble-like flow ( [Lee et al.|2000| ), where the incli¬ 
nation then enters as a free parameter. It is defined as the angle 
between the outflow direction and the line of sight ( jCabrit &] 


|Bertout|[l990| /=0° is pole on). Small radial velocities are ex¬ 


pected for outflows which lie in the plane of the sky. Therefore a 
correction factor for inclination Cf is applied in the calculations. 
In Table [T] the correction factors from [Downes & Cabrit| ( |2007] ) 
are tabulated; these correction factors come from detailed out¬ 
flow modeling and synthetic observations of the model results. 
Moreover, we note that these correction factors include correc¬ 


tion for missing mass within + 2 km s“^ from the source veloc¬ 
ity. The correction factors have been applied to the outflow rate, 
force and luminosity as listed in Tables and The velocity, 
as a measured parameter, is not corrected for inclination. The 
inclination angles are estimated from the outflow maps as fol¬ 
lows: if the outflow lobes are overlapping, the outflow is likely 
very inclined. If the outflow shows low-velocity line wings but a 
large extent on the sky, the inclination is very likely low. In this 
way each outflow is classified individually, and divided into in¬ 
clination bins at 10°, 30°, 50°, and 70°. Our estimates are listed 


available ^abrit & Bertout 

1992[ [Gueth et al. 


Bourke 

et al.|1997[ Hogerheijde et a] 

L|1997 [Micono et sl. 


[Brown 


& Chandler[[1999t[Lommen et al. [[20081 [Tobin et al. [[2668 [ [van 


Kempen et al.|2009b| ), e xcept for IRAS 15398 for w hich we And 

a larger inclination than [van Kempen et al.[ ( [2009^ . Our inclina¬ 
tion of IRAS 15398 is consistent with newer values from (Oya 
[et al.[2014| ). Although the method for determining the outflow in¬ 
clinations is subjective, the inclinations agree with literature val¬ 
ues where available, which lends some credibility to the method, 
and we estimate that the uncertainty is 30°. That is, the correc¬ 
tion introduces a potential systematic error of up to a factor of 2 
in the outflow parameters. 

The resulting maps of all sources are presented in Figs. 
and for 6-5 and 3-2, respectively, where blue and red 
contours show the blue- and red-shifted outflow lobes, respec- 
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Table 1. Inclination correction factors. 


in 

10 

30 

50 

70 

Oi 

1.2 

2.8 

4.4 

7.1 


Notes. Line-of-sight i nclinations, where i-0° indicates pole-on 
( [Downes & Cabrit|2007| ). 


lively. The velocity limits are summarized in Table [A.lj in the 
Appendix. A few maps cover only the central ~2'x2', specifi¬ 
cally the three Class 0 sources NGC 1333-IRAS 2A, L723mm, 
L1527, and the two Class I sources Elias 29 and L1551-IRS5. 
Source-by-source outflow and intensity maps obtained from the 
CO 6-5 and 3-2 data are presented in Figs. |A.2 
line Appendix. 

2.4. CO maps 

The '^CO 6-5 and 3-2 transitions were mapped around 
the central ~TxT region, corresponding to typically 
-10^ AU X 10"^ AU. The total integrated intensity i s mea¬ 
sured for all the sources and presented in Table C.1-26 of jYildizj 
jet al.j (|2013|) for the source positions. All maps are presented 

and |A.6| and as spectral maps in 
Figs.|A.8|A.13|in the Appendix. 


as contour maps in Figs. A.5 


A.4 in the On- 


3. Results 


3.1. Outflow morphology 

All sources show strong outfiow activity in both CO transitions, 
7=6-5 and 3-2, as is evident from both the maps and spectra 
(Figs. 1313 and Figs. |A.lf[A.4| ). The advantage of the CO 6-5 
maps is that they have higher spatial resolution by a factor of 2 
than the CO 3-2 maps. On the other hand, the CO 3-2 maps have 
the advantage of higher S /N than the CO 6-5 maps by typically 
a factor of 4 in main beam temperature. 

Most sources show a clear blue-red bipolar structure. In a few 
cases only one lobe is observed. Specific examples are TMCIA, 
which shows no red-shifted outfiow lobe, and HH 46, which has 
only a very small blue-shifted outfiow lobe. One explanation is 
that these sources are at the ed ge of the cloud and that th ere is 
no cloud material to run into ( van Kempen et al]|20Q9b| ). For 
F723mm, NGC 1333-IRAS 2A and BHR71, two outflows are 
driven by two independent pr otostars ( [Lee et aL]|2002t jParise 
jet al.|2Q06l|Codella et al.|2Q14| ) and both outflows are detected in 
our CO 3-2 maps. In CO 6-5, only one outflow shows up toward 
F723mm and NGC 1333-IRAS 2A, whereas both outflows are 
seen toward BHR71. 

Visual inspection shows that the Class 0 outflows are more 
collimated than their Class I counterparts as expected (e.g., |Arce 
et al.|2007| ). The length of the outflows can be quantified for most 


of the sources. Rqo is defined as the total outflow extent assum¬ 
ing that the outflows are fully covered in the map. Rqo is mea¬ 
sured separately for the blue and red outflow lobes as the pro¬ 
jected size, with sometimes significantly different values. Rqo as 
measured from CO 6-5 is applied to CO 3-2 in the cases where 
the CO 6-5 maps are larger than their 3-2 counterparts. To¬ 
ward some sources, e.g., DK Cha and NGC 1333-IRAS 4B, the 
blue and red outflow lobes overlap, likely because the outflows 
are observed nearly pole on. In other cases the outflow lobes 
cannot be properly isolated from nearby neighboring outflow 
lobes. Such a confusion is most pronounced in Ophiuchus (e.g.. 
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Fig. 4. Histogram of total Rqo (blue- and red-shifted outflows com¬ 
bined) is shown for Class 0 (red) and Class I (blue) sources. (Rqo is 
not corrected for inclination.) 



Fig. 5. Rqo is plotted against Riok. the radius of the modeled envelope 
within 10 K radius. The black line is for Rqo = ^iok, showing that al¬ 
most all sources follow Rqo > Riok and that Rqo is larger for Class 0 
than Class I sources. 


GSS30-IRS1). In those cases, Rqo could not be properly esti¬ 
mated and the estimated value is a lower limit. Figure]^ shows a 
histogram of total Rqo for Class 0 and I sources. Class 0 sources 
show a nearly flat distribution across the measured range of ex¬ 
tents, whereas few Class I sources show large outflows (L1551 
is a notable exception). In Fig.[^ Rqo is plotted against Riok, the 
radius of the modeled envelope within a 10 K radius. The out¬ 
flowing gas typically extends to much greater distances than the 
surrounding envelope and thus influences the surrounding cloud 
material directly. 


3.2. Outflow parameters 


In the following, different outflow parameters, including mass, 
force and luminosity, are measured. These parameters have pre¬ 
viously been determined from lower-/ lines for several young 


stellar objects (e.g., jCabrit & Bertout 1992[ Bontemps et al. 

1996 

Hogerheijde et al.|1998 [Hatchell et al.|2007 

Curtis et al. 

2010 

van der Marel et al.|2013l Dunham et al.|20] 

and more 

recently from CO 6-5 by [van Kempen et al. (2009 

3) andjYildiz 


results are listed in Tables |2] and |3l Uncertainties in the methods 
are discussed extensively in [van der Marel et ar] ( |2013| ). 
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Table 2. Outflow properties of the red and blue outflow lobes of Class 0 sources. 


Source 

Trans. 

Inclination 

n 

Lobe 

Rco‘‘ 

[AU] 

f a,b 

^dyn 

[103 yr] 

outflow 

[Mo] 

M‘‘-‘ 

[Moyr-'l 

Fco‘‘^ 

[Mq yr ‘km s *] 

ikin'*’** 

[LgI 

L1448MM 

CO 3-2 

50 

Blue 

5.9x10'^ 

5.5 

9.0x10-3 

7.2x10-3 

2.0x10-3 

2.8x10*3 




Red 

5.9x10'* 

9.7 

6.2x10-3 

2.8x10-3 

1.7x10-3 

2.3x10** 

NGC1333-IRAS 2A 

CO 6-5 

70 

Blue 

1.4x10'* 

2.9 

7.9x10-3 

2.0x10-3 

3.4x10-“* 

1.7x10-* 




Red 

1.4x10'* 

3.9 

2.2x10-3 

4.0x10-3 

2.0x10-3 

1.7x10** 


CO 3-2 

70 

Blue 

2.4x10'* 

4.8 

8.5x10-3 

1.3x10-“* 

2.6x10-3 

1.2x10** 




Red 

2.4x10“* 

6.4 

6.9x10-3 

7.7x10-3 

4.8x10-3 

5.4x10** 

NGC1333-IRAS 4A 

CO 6-5 

50 

Blue 

2.5x10“* 

5.3 

8.1x10-3 

6.7x10-3 

1.5x10-“* 

5.6x10-3 




Red 

3.5x10“* 

8.4 

1.9x10-3 

9.9x10-3 

5.4x10-“* 

5.4x10-* 


CO 3-2 

50 

Blue 

2.8x10“* 

6.1 

2.1x10-3 

1.5x10-3 

3.5x10-“* 

1.6x10-* 




Red 

3.9x10“* 

9.3 

2.5x10-3 

1.2x10-3 

1.7x10-3 

1.8x10** 

NGC1333-IRAS 4B 

CO 6-5 

10 

Blue 

2.4x10^ 

0.6 

8.2x10-“* 

1.6x10-3 

3.2x10-3 

7.4x10-3 




Red 

1.2x103 

0.4 

7.3x10-“* 

2.2x10-3 

1.6x10-“* 

1.8x10-* 


CO 3-2 

10 

Blue 

3.5x103 

0.8 

8.3x10-“* 

1.3x10-3 

2.9x10-3 

1.1x10-3 




Red 

2.4x103 

0.9 

2.7x10-3 

3.6x10-3 

1.9x10-“* 

1.8x10-* 

L1527 

CO 6-5 

70 

Blue 

1.5x10“* 

9.1 

2.3x10-3 

1.8x10-3 

3.1x10-3 

9.9x10-3 




Red 

1.1x10“* 

6.5 

2.5x10-3 

2.7x10-3 

1.1x10-“* 

7.5x10-3 


CO 3-2 

70 

Blue 

3.2x10“* 

20.6 

1.0x10-3 

3.5x10-3 

6.1x10-3 

2.0x10-3 




Red 

1.1x10“* 

6.5 

9.0x10-3 

9.8x10-3 

3.8x10-“* 

2.6x10-* 

CedllO-IRS4 

CO 6-5 

30 

Blue 

3.8x103 

4.2 

2.7x10-“* 

1.8x10-' 

2.1x10-3 

4.6x10-“* 




Red 

3.8x103 

4.7 

2.5x10-“* 

1.5x10-3 

4.1x10-3 

2.0x10-3 

BHR71 

CO 6-5 

70 

Blue 

4.4x10“* 

13.4 

3.4x10-3 

1.8x10-3 

7.7x10-“* 

6.2x10-* 




Red 

4.0x10“* 

8.5 

6.9x10-3 

5.8x10-3 

7.7x10-“* 

3.3x10-* 

IRAS 15398 

CO 6-5 

30 

Blue 

2.6x103 

1.4 

3.4x10-“* 

1.7x10-3 

6.4x10-3 

1.4x10-3 




Red 

2.6x103 

1.2 

2.7x10-“* 

1.5x10-3 

2.6x10-3 

2.0x10-3 


CO 3-2 

30 

Blue 

3.2x103 

1.8 

4.4x10-“* 

1.8x10-3 

9.2x10-3 

2.5x10-3 




Red 

2.0x103 

0.9 

2.5x10-“* 

1.9x10-3 

2.8x10-3 

2.0x10-3 

L483MM 

CO 6-5 

70 

Blue 

1.2x10“* 

5.2 

4.2x10-3 

5.7x10-3 

6.7x10-3 

1.6x10-3 




Red 

1.0x10“* 

4.4 

3.4x10-3 

5.4x10-3 

2.1x10-“* 

1.5x10-* 


CO 3-2 

70 

Blue 

1.4x10“* 

6.2 

7.0x10-3 

8.0x10-3 

7.7x10-3 

1.6x10-3 




Red 

1.0x10“* 

4.4 

8.5x10-3 

1.4x10-3 

5.1x10-“* 

3.4x10-* 

Ser-SMMl 

CO 6-5 

50 

Blue 

3.4x10“* 

8.4 

1.6x10-3 

8.2x10-3 

1.5x10-“* 

5.7x10-3 




Red 

1.8x10“* 

3.9 

1.2x10-3 

1.4x10-3 

8.7x10-“* 

9.8x10-* 


CO 3-2 

50 

Blue 

3.4x10“* 

8.6 

6.4x10-3 

3.3x10-3 

6.7x10-“* 

2.8x10-* 




Red 

1.8x10“* 

3.9 

3.3x10-3 

3.7x10-3 

2.3x10-3 

2.7x10** 

Ser-SMM4 

CO 6-5 

30 

Blue 

1.8x10“* 

4.6 

2.4x10-3 

1.5x10-3 

2.5x10-“* 

9.3x10-3 




Red 

1.8x10“* 

7.3 

2.8x10-3 

1.1x10-3 

5.9x10-“* 

5.8x10-* 


CO 3-2 

30 

Blue 

1.8x10“* 

4.6 

1.6x10-* 

9.9x10-3 

2.0x10-3 

8.3x10-* 




Red 

1.8x10“* 

7.6 

1.3x10-* 

4.7x10-3 

2.8x10-3 

2.9x10** 

Ser-SMM3 

CO 6-5 

50 

Blue 

4.6x103 

1.0 

6.9x10-3 

3.1x10-3 

6.0x10-“* 

2.6x10-* 




Red 

4.6x103 

1.6 

3.0x10-3 

8.1x10-3 

4.9x10-“* 

5.4x10-* 


CO 3-2 

50 

Blue 

4.6x103 

1.0 

2.7x10-3 

1.2x10-“* 

2.4x10-3 

1.0x10** 




Red 

4.6x103 

1.6 

1.1x10-3 

3.0x10-3 

1.8x10-3 

1.9x10** 

B335 

CO 6-5 

70 

Blue 

6.2x103 

3.4 

4.7x10-“* 

9.9x10-’ 

2.3x10-3 

9.4x10-3 




Red 

8.8x103 

4.8 

1.3x10-3 

1.9x10-3 

9.0x10-3 

7.6x10-3 


CO 3-2 

70 

Blue 

1.0x10“* 

5.3 

3.7x10-3 

4.9x10-3 

1.3x10-“* 

6.2x10-3 




Red 

7.5x103 

4.1 

5.4x10-3 

9.3x10-3 

4.7x10-“* 

4.3x10-* 

L723MM 

CO 6-5 

50 

Blue 

1.2x10“* 

4.1 

6.0x10-3 

6.6x10-3 

1.8x10-“* 

1.0x10-* 




Red 

1.2x10“* 

3.8 

7.5x10-3 

8.5x10-3 

5.5x10-“* 

6.3x10-* 


CO 3-2 

50 

Blue 

1.8x10“* 

6.0 

3.0x10-3 

2.2x10-3 

6.5x10-“* 

3.7x10-* 




Red 

1.8x10“* 

5.8 

4.2x10-3 

3.2x10-3 

2.2x10-3 

2.8x10** 

LI 157 

CO 3-2 

70 

Blue 

4.4x10“* 

16.8 

1.2x10-* 

4.9x10-3 

5.0x10-“* 

1.9x10-* 




Red 

5.2x10“* 

14.1 

1.5x10-* 

7.3x10-3 

3.2x10-3 

3.1x10** 


Notes. Outflow extents and outflow masses are not corrected for inclination. Dynamical times cale . Constant temperature of 75 K is 
assumed for both CO 6-5 and CO 3-2 calculations. Corrected for inclination as explained in Sect. 3.2 Mass outflow rate ^ Outflow force 
Kinetic luminosity. 


3.2.1. Outflow mass 

One of the most basic outflow parameters is the mass. The in¬ 
ferred mass depends on three assumptions: the line opacity, the 


distribution of level populations, and the CO abundance with re¬ 
spect to H 2 . In the following, we assume that the line wings are 
optically thin, as has been demonstrated observationally for CO 
6-5 for a few sources with massive outflows (e.g., NGC 1333- 
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Table 3. Outflow properties of the red and blue outflow lobes of Class I sources. 


Source 

Trans. 

Inclination 

n 

Lobe 

Rco^ 

[AU] 

f a,b 

^dyn 

[103 yr] 

outflow 

[Mo] 

M‘‘-‘ 

[Mo yr-‘] 

Fco‘‘^ 

[M© yr ^km s 

[Lo] 

L1489 

CO 6-5 

50 

Blue 

3.5x10^ 

1.2 

4.0x10-3 

1.5x10-' 

1.8x10-3 

3.2x10-“ 




Red 

2.1x103 

1.3 

1.3x10-“ 

4.6x10-3 

2.3x10-3 

2.0x10-3 


CO 3-2 

50 

Blue 

3.5x103 

1.2 

6.9x10-“ 

2.5x10-'' 

3.9x10-3 

1.3x10-3 




Red 

2.1x103 

1.3 

7.1x10-“ 

2.5x10-'' 

1.2x10-“ 

1.1x10-' 

L1551-IRS5 

CO 3-2 

70 

Blue 

l.TxlO-* 

8.2 

7.4x10-3 

6.4x10-'’ 

9.3x10-3 

2.7x10-3 




Red 

l.TxlO-* 

6.7 

9.6x10-3 

1.0x10-3 

4.2x10-“ 

3.1xl0-‘ 

TMRl 

CO 6-5 

50 

Blue 

4.9x103 

2.9 

1.2x10-“ 

1.8x10-3 

2.0x10-3 

4.9x10-“ 




Red 

3.5x103 

4.5 

2.4x10-“ 

2.3x10-3 

8.0x10-3 

4.8x10-3 


CO 3-2 

50 

Blue 

4.9x103 

3.0 

2.6x10-“ 

3.8x10-3 

5.1x10-3 

1.4x10-3 




Red 

3.5x103 

4.5 

5.8x10-“ 

5.7x10-3 

2.0x10-3 

1.3x10-3 

TMCIA 

CO 6-5 

50 

Blue 

5.6x103 

1.4 

2.3x10-“ 

7.2x10-3 

1.3x10-3 

7.6x10-3 




Red 

2.1x103 

1.8 

4.0x10-'’ 

9.5x10-5 

3.7x10-3 

2.4x10-“ 


CO 3-2 

50 

Blue 

5.6x103 

1.6 

2.8x10-3 

7.8x10-'' 

1.1x10-“ 

3.7x10-3 




Red 

1.7x103 

1.5 

1.4x10-“ 

4.1x10-3 

1.8x10-3 

1.4x10-3 

TMCl 

CO 6-5 

50 

Blue 

3.5x103 

1.2 

1.3x10-“ 

4.7x10-3 

3.2x10-3 

4.0x10-“ 




Red 

4.9x103 

1.6 

3.8x10-“ 

1.1x10-'' 

5.0x10-3 

4.4x10-3 


CO 3-2 

50 

Blue 

3.5x103 

1.2 

5.6x10-“ 

2.0x10-'' 

2.7x10-3 

7.9x10-3 




Red 

2.1x103 

0.7 

1.4x10-3 

8.9x10-'' 

4.2x10-“ 

3.7xl0-‘ 

HH46-IRS 

CO 6-5 

50 

Blue 

1.1x10“ 

9.9 

2.6x10-3 

1.2x10-'’ 

1.2x10-3 

2.9x10-3 




Red 

2.5x10“ 

7.9 

3.2x10-3 

1.8x10-3 

7.7x10-“ 

6.5xl0-‘ 


CO 3-2 

50 

Blue 

1.6x10“ 

13.6 

2.2x10-3 

7.2x10-3 

2.5x10-“ 

1.7xl0-‘ 




Red 

2.5x10“ 

7.9 

2.2x10-3 

1.2x10-3 

8.1x10-“ 

9.3xl0-‘ 

DK Cha 

CO 6-5 

10 

Blue 

1.8x103 

1.6 

1.8x10-“ 

1.3x10-3 

6.6x10-3 

7.3x10-3 




Red 

1.8x103 

0.9 

1.1x10-“ 

1.4x10-3 

2.5x10-3 

4.4x10-“ 

GSS30-IRS1 

CO 6-5 

30 

Blue 

1.5x10“ 

5.5 

1.5x10-3 

7.9x10-3 

5.1x10-3 

1.1x10-3 




Red 

1.5x10“ 

4.9 

8.9x10-3 

5.1x10-3 

2.0x10-“ 

1.6xl0-‘ 


CO 3-2 

30 

Blue 

1.5x10“ 

5.5 

2.1x10-3 

1.1x10-3 

6.0x10-3 

8.8x10-3 




Red 

1.5x10“ 

4.9 

2.4x10-3 

1.4x10-3 

4.6x10-“ 

3.0xl0-‘ 

Elias 29 

CO 6-5 

30 

Blue 

7.5x103 

3.1 

6.4x10-“ 

5.7x10-3 

4.4x10-3 

1.2x10-3 




Red 

5.0x103 

1.7 

6.3x10-“ 

1.0x10-3 

3.9x10-3 

2.7x10-3 


CO 3-2 

30 

Blue 

7.5x103 

3.6 

1.4x10-3 

1.1x10-3 

6.6x10-3 

1.0x10-3 




Red 

7.5x103 

3.3 

1.8x10-3 

1.5x10-3 

5.7x10-3 

3.9x10-3 

Oph-IRS63 

CO 6-5 

50 

Blue 

3.8x103 

1.6 

1.0x10-“ 

2.8x10-3 

4.3x10-3 

2.1x10-3 




Red 

3.8x103 

4.2 

8.6x10-3 

9.0x10-* 

2.1x10-3 

8.3x10-“ 


CO 3-2 

50 

Blue 

8.8x103 

3.7 

7.0x10-“ 

8.4x10-3 

5.6x10-3 

1.5x10-3 




Red 

5.0x103 

7.4 

5.0x10-“ 

3.0x10-3 

5.2x10-3 

2.2x10-3 

RN091 

CO 6-5 

50 

Blue 

3.1x103 

1.0 

2.5x10-“ 

1.1x10-3 

2.9x10-3 

2.1x10-3 




Red 

1.9x103 

2.5 

7.3x10-3 

1.3x10-3 

1.1x10-3 

1.8x10-“ 


CO 3-2 

50 

Blue 

6.2x103 

2.0 

3.5x10-3 

7.7x10-3 

1.0x10-“ 

4.1x10-3 




Red 

1.9x103 

2.5 

2.3x10-“ 

4.0x10-3 

2.3x10-3 

2.6x10-“ 


Notes. Outflow extents and outflow masses are not corrected for inclination. Dynamical time scale . Constant temperature of 75 K is 
assumed for both CO 6-5 and CO 3-2 calculations. Corrected for inclination as explained in Sect. 3.2 Mass outflow rate ^ Outflow force 
Kinetic luminosity. 


IRAS 4A, [Yildiz et al.||2012| ). CO 3-2 emission is also as¬ 
sumed optically thin in the following, although that assumption 
may not be fully valid (see discussion below). The level popu¬ 
lations are assumed to follow a Boltzmann distribution with a 
single temperature, Tgx. F inally, the ab undance ratio is taken as 
[H2/'^CO]=1.2 x 104, as in |Yildiz et all ( |2012| . 

The upper level column density per statistical weight in a 
single pixel (4'.'5x4'.'5 for CO 6-5, 7'.'5x7'.'5 for CO 3-2) is cal¬ 
culated as 


/3v^ f r^bdV 

gu ^ul gu 


The total CO column density in a pixel, A^totaU is 

Votai = ; ( 2 ) 

Q(T) is the partition function corresponding to a specific exci¬ 
tation temperature, Tgx, which is assumed to be 75 K for both 
CO 3-2 and CO 6-5 ob servations ( [van Kempen et al.||2009bt 
Yildiz et ar]|2Q12| |2Q13| ). Changing Tex by +30 K changes the 
inferred column densities by only 10-20%. 

The mass is calculated as 


( 1 ) 


outflow ~ /^H 2 ^ 


The constant p is /hc^ = 1931 cm~^ (GHz^ K km)“^ The re¬ 
maining parameters are for the specific transition, where v is the 
frequency, Aui is the Einstein A coefficient and g^=2J-\-l. 
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H2 


12CO 


^ A/total,; 


(3) 


where the factor includes the contribution of helium 

( Kauffmann et al.|[2008| ) and mn is the mass of the hydrogen 
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atom. A is the surface area of one pixel j. The sum is over all 
outflow pixels. 

The mass may be underestimated if the line emission is 
optically thick. data exist toward most outflows (see above) 
but the S /N of these data is typically too low to properly mea¬ 
sure the opacity in the line wings, except for at the source po¬ 
sition where the signal naturally is the strongest. The line 
wings do not extend beyond the inner velocity limits. NGC1333- 
IRAS 4A is one of the few sources where line w ings are de¬ 
tected in at the outflow positions (Fig. 11 in 


Yildiz et al. 


2012), and it is clear that at the velocity ranges considered here, 
the line emission is optically thin (r < 1); the same is true for 
the outflows studied by [van der Marel et aL| ( |2013| ) in CO 3-2 
emission in Ophiuchus (their Fig. 4), where deep pointed ob¬ 
servations of 3-2 were required to measure the opacity. 
That study concluded that the opacity does not play a signifi¬ 
cant role when determining the outflow parameters. Similarly, 
[Dunham et al.| p014| ) conclude that CO 3-2 may be optically 
thick at velocities less than 2 km s“^ offset from the source ve¬ 
locity, velocities which are excluded from our analysis because 
of the risk of cloud contamination. Potentially more problem¬ 
atic is the missing mass at low velocities. The missing mass is 
moving close to the systemic velocity and it is not possible to 
disentangle this mass from the surrounding cloud material, an 
effect which may introduce a typical uncertainty of a factor of 
2-3 ( [Downes & Cabrit||2007| ). However, the correction factors 
derived by the same authors and implemented here account for 
that missing mass. 6-5 emission will be less affected by 
this confusion than the 3-2 emission, simply because of 
the different excitation conditions required. 


3.2.2. Outflow force 

One of the most important outflow parameters is the outflow 
force, Fco- The best method for computing the outflow force 
is still debated and the results suffer from ill-constrained obser¬ 
vational parameters, such as inclination, i. [van der Marel et aL 
( 2013[ ) compare seven different methods proposed in the litera¬ 
ture to calculate outflow forces. The “separation method” (see 
below) in their paper is found to be the preferred method, which 
is less affected by the observational biases. The method can also 
be applied to low spatial resolution observations or incomplete 
maps. Uncertainties are estimated to be a factor of 2-3. 

In the following, the outflow force is calculated separately 
for the blue- and red-shifted lobes, only including emission 
above the 3cr level. The mass is calculated for each channel sep¬ 
arately and multiplied by the central velocity of that particular 
channel. The integral runs over velocities from Vin to Vout- They 
are then summed and the sum is over all pixels j in the map with 
outflow emission. The outflow force is calculated for the red- and 
blue-shifted outflow lobes separately. This method is formulated 
as: 

VmaxZJ/WOV'dy'j. 

Fco = Ci ---^, (4) 

^CO 

where c/ is the inclination correction (Table [T]), and Rqo is the 
projected size of the red- or blue-shifted outflow lobe. The out¬ 
flow force is computed separately from the CO 3-2 and 6-5 
maps of the same source (see Tables and |^. 

The difference in outflow force between the red and blue out¬ 
flow lobes ranges from ~1 up to a factor of 10. For sources with 
a low outflow force such as Oph IRS63 (<10“^ Mq yr“^ km s“^) 
this is a result of differences in the inferred outflow mass per 





Fig. 6. Outflow forces {left) and outflow masses {right), calculated from 
CO 6-5 and 3-2 emission are compared for Class 0 and I sources. Green 
lines are for a ratio of 1. 



Fig. 7. Histograms of calculated total outflow force Fco are shown for 
Class 0 (red) and Class I (blue) sources. 


lobe, which, in these specific cases, is primarily a result of low 
S IN. In these cases, the overall uncertainty on the outflow force 
is high, up to a factor of 10. In other cases, such as HH46 as 
mentioned above, there is a real asymmetry between the differ¬ 
ent lobes which is caused by a difference in the surrounding en¬ 
vironment. In the following, only the sum of the outflow forces 
of both lobes as measured from each outflow lobe will be used. 

Figure shows how the outflow forces and outflow masses 
calculated from CO 3-2 and 6-5 differ. For strong outflows, 
there is a factor of a few difference in the two calculations, with 
differences up to an order of magnitude for the weaker outflow 
sources. Although the CO 6-5 emission suffers less from opac¬ 
ity effects and so recovers more emission / mass at lower veloci¬ 
ties, this effect is overwhelmed by the lower S /N of the CO 6-5 
emission. The fact that the masses and outflow forces derived 
from the 6-5 data are systematically lower than those from the 
3-2 d ata is likely due to the same effect ( [van der Marel et aL 
2013[). Moreover, if CO 6-5 traces slightly warmer gas than CO 
3-2 ( Yildiz et al.|2Q13| ) then the mass traced by this line will be 
lower than that traced by CO 3-2. Both effects work to system¬ 
atically lower the CO 6-5 masses, which in turn leads to lower 
outflow forces. 

Figure [^displays Fco from CO 6-5 for Class 0 and Class I 
sources separately. Generally, Class 0 sources have higher out¬ 
flow forces and are thus more powerful than their Class I coun¬ 
terparts ( [Bontemps et al.|1996| ). The Class I source with an ex¬ 
ceptionally high outflow force is HH46. 
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Table 4. Median values of the outflow parameters. 



^outflow 

[Mo] 

M 

[Moyr-*] 

Leo 

[Mokm^ yr"^] 

^kin 

[iol 

CO 6-5 

Class 0 

9.8x10-^ 

1.5x10-3 

6.9x10-* 

6.0x10-* 

Class I 

3.4x10-* 

8.4x10-*^ 

2.8x10-3 

1.2x10-* 

Total 

2.2x10-^ 

1.0x10-3 

1.4x10-* 

1.9x10-* 

CO 3-2 

Class 0 

7.2x10-^ 

5.4x10-3 

2.9x10-3 

2.9x10** 

Class I 

3.1x10-3 

3.0x10-3 

1.4x10-* 

1.4x10** 

Total 

1.7x10-3 

3.3x10-3 

5.2x10-* 

1.9x10** 


3.3. Other outflow parameters 

Other outflow parameters that characterize the outflow activity 
are the dynamical age, ^dyn, mass outflow rate, Moutfiow. and ki¬ 
netic luminosity, L^in- 

Assuming that the outflow moves with a constant velocity 
over the extent of the outflow, the dynamical age is determined 
as 


^dyn ~ 


krnax 


(5) 


This age is a lower limit on the age of the protostar ( [Curtis et aL] 


2010) if the outflowing material is decelerated, e.g., through in¬ 


teractions with the ambient surrounding material. On the other 
hand, the outflow may be signiflcantly younger since the veloc¬ 
ities of the central jet that drives the molecular outflow are typi¬ 
cally higher than 100 km s“^ and what is observed in these colder 
low-/ CO lines may just be the outer shell which is currently 
undergoing acceleration, not deceleration. See, e.g., [Downes &| 
[Cabrit[ ( [2()07[ ) for a more complete discussion. The outflow mass 
loss rate is computed according to 


)l‘/f ^outflow 

^outflow ~ ~ • (6) 

^dyn 

The kinetic luminosity is given by 

^kin — 2'^Col^max • (7) 

Outflow parameters of Fqo, and Lkin with inclination cor¬ 
rections are presented in Tables [^and|^ However, Moutflow. 

^dyn, and Vmax are not corrected for inclination, since they are 
measured quantities. The median values of the results are given 
in Table [U 


3.4. Correlations 


Most previous studies of the outflow force were done using CO 


1-0, 2-1, or 3-2 (e.g., Cabrit & Bertout[[1992 

Bontemps et al. 

19961 Hogerheijde et al. [1998 

Hatchell et al. 

[20071 van Kem- 

pen et al.[[2009c[ Dunham et a 

l.[[2014[). The opacity decreases 


with excitation, as suggested by, e.g., the observations reported 
Dunham et al.[ ( [2014[ ), but without targeted, deep surveys of 


m 

13 


CO, it is difficult to quantify how much the CO column den¬ 
sity is underestimated. Furthermore, cloud or envelope emission 
may contribute to the emission at the lowest outflow velocities 
at which the bulk of the mass is flowing. With our CO 6-5 ob¬ 
servations, some of the above-mentioned issues can be avoided, 
or their effects can be lessened. Thus, it is important to revisit 
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the correlations of outflow force with bolometric luminosity and 
envelope mass using these new measurements. 

In Fig. [8] Fco is plotted against Lboi, Mgnv, and Moutflow, 
where the Fqo and Moutflow values are taken from the CO 6-5 
data. The best fit between Fqo and Lboi is shown with the green 
line corresponding to 

log(Fco) = -(4.71 + 0.02) + (1.13 + 0.37)log(Lboi). (8) 

Outflows from Class 0 and Class I sources are well-separated; 
Class 0 sources show more powerful outflows compared to 
Class I sources of similar luminosity. The Pearson correlation 
coefficients are r=0.62, 0.83, and 0.64 for all sources. Class 0, 
and Class I sources, corresponding to confidences of 2.9, 2.9, 
and 1.9(T, respectively. 

The best fit between Fqo and Menv is described as 


log(Fco) = -(3.95 + 0.37) + (1.24 + 0.21) log(Menv) (9) 

and Pearson correlation coefficients are r=0.81, 0.82, and 0.56 
(3.8, 2.8 and l.lcr) for all sources. Class 0, and Class I, re¬ 
spectively. Since early Class 0 sources have higher accretion 
rates their outflow force is much higher than for the Class I 
sources (see, e.g., [Bontemps et ffi^[1996[ for a full discussion). 
Finally, as expected, a strong correlation is found between Fqo 
and Moutflow with a Pearson correlation coefficient of r=0.92 for 
all sources (4.3cr), not surprisingly since Fqo is nearly propor¬ 
tional to Moutflow- The best fit is described as 


log(Fco) = -(1.71 ± 0.02) -r (0.88 + 0.62) log(Mco). 


( 10 ) 


Previously, [Bontemps et al.[ ( [1996[ ) surveyed 45 sources us¬ 
ing CO 2-1 observations with small-scale maps. In Fig. the 
blue and green dashed lines of Fqo vs. Lboi and Mgnv show the 
fit results from their Figs. 5 and 6 ( [Bontemps et al.|1996] ). Since 
their number of Class I sources is higher than Class 0 sources, 
the fit was only done for Class I sources in Leo vs. Lboi. In Fig.[^ 
the blue solid line only shows the fit for Class I sources and the 
correlation is described by. 


log(Lco) = -(5.14 + 0.29) + (0.98 + 0.55) log(Lboi). 


( 11 ) 


In the Leo vs. Mgnv plot, the fits are shown as green lines for the 
entire sample. The [Bontemps et al.[ ( [1996[ ) sample is weighted 
toward lower luminosities (< 10 Lboi), where our Leo measure¬ 
ments from the CO 6-5 data follow their relation for Class I 
sources obtained from 2-1 data, but with a shift to a factor of 
a few higher values of Leo- However, given the scatter in the 
results for low Lboi sources, this difference is hardly significant. 

Examining the same outflow parameters measured using the 
CO 3-2 transition, and their correlation with the same outflow 
parameters, a similar picture arises (Fig. [A.7[ ). However, for the 
sources in our sample, the correlations follow the same trend but 
they are somewhat weaker. In particular, the correlation with Lboi 
is at the ~2.7(T level, whereas the correlation with Mgnv is S.lcr. 
Although the measured values of, e.g., Leo, hh out the same 
parameter space as when the measurements are done with CO 
6-5, the scatter is larger. The scatter remains on the order of one 
order of magnitude, which is similar to the scatter reported in the 
literature (e.g., [Bontemps et al.|1996[ ), but because of the limited 
source sample (20 sources with Leo measurements) it is difficult 
to compare these 3-2 measurements with what is presented in 
the literature. 
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Fig. 8. Correlations of Fqo with Lboi, Mgnv, and Moutflow. where Fco is determined from the CO 6-5 data. Red and blue symbols indicate Class 0 
and Class I sources, respectively. The g reen solid line is the fit to all values and the blue solid line is the fit to the Class I sources alone. Blue and 
green dashed lines are the best fits from |Bontemps et ar] ( |1996| ). 


3.5. Radiative feedback from UV heating 

The quiescent gas is traced by the narrow (FWHM <1 km s“^) 
6-5 emission, which has been mapped over a 1' region 
aro und th e source position. As the contour maps in Figs. |A.5| 
and |A.6| show, the emission is strongly centrally concentrated 
and does not extend beyond the mapped region except for spe¬ 
cial cases like NGC 1333 IRAS 4A ( |Yildiz et al.||MT2| ). The 
observed emisison has two contributions: (/) the dense envelope 
heated ‘passively’ by the luminosity of the protostar, i.e., the 
dust in the envelope absorbs all the protostellar luminosity and 
is heated by it, and this temperature is then transferred to the 
gas through gas-dust collisions; (ii) the gas heated by UV pho¬ 
tons created by protostellar accretion or by shocks in the outflow, 
and escaping from the immediate protostellar surroundings, for 
example through outflow cavities, to larger distances. Here the 
temperature of the gas is higher than that of the dust. 

The first component has been modeled by [Kristensen et al. 
( |2012| ) for all our sources and dust temperatures in excess 
of 10 K are typically found out to between 2.5x10^ up to 
1.5x10'^ AU from the sources. There is evidence that the dust 
may be further heated on large sc ales by the UV photons gen- 
erated by the accreting protostar ( [Hatchell et al.||2Qf3{ [Sicilia- 
I Aguilar et al.||2Q13] ). We here quantify the second component, 
which is the gas with temperatures higher than that of the dust, 
in excess of the passively heated envelope. This second mecha¬ 
nism operates on larger scales and is most rel evant to preventing 
further collapse or fragmentation of the core dOffner et al. 120091 

[ 2 ^ . -- 

To isolate this second component, the method outlined in 


Yildiz et al. 


( 2012| is used. The 6-5 envelope emission 


(component ( 0) is modeled using the temperature and density 
profiles from Kristensen et al. (2012) together with the 
constant abundance results provided in | Yildiz etH] p013[ Ta¬ 
ble 5). For the three NGC 1333 sources, drop abundance profiles 
are used in which CO is frozen out in some part of the envelope; 
for NGC 1333-IRAS 2A the results from [Yildiz eTak] dMIO] ) 
are taken, whereas for NG C 1333-IRAS 4A an d NGC 1333- 
IRAS 4B the models from | Yildiz et~ar] ( |2012| ) are adopted. 
These abundances are then multiplied by t he [^^C]/[^^0] 
abundance ratio of 8.5 ( Langer & Penzias||l990| ) and the 
emission is computed using the non-LTE excitation line radia¬ 
tive transfer code RATRAN ( [Hogerheijde & van der Tak|2000| ). 
The turbulent width for all the model spectra is taken as 


0.9 km s-\ except for NGC 13 33-IRAS 4A and N GC 1333- 
IRAS 4B where the values from [Yildiz et~ar] ( |2Q12| ) are used. 
The resulting emission map is convolved with the relevant ob¬ 
serving beam. 

Figures A.8|A.13| present 7x7 pixel maps (~30" x 30"; 1 
pixel = 4.5") around the central position of each source in the 
and 6-5 transitions. The modeled envelope emission 
(component (/)) is shown as red lines overplotted on the 
maps. The right-most panels present the difference between the 
model envelope emission and the observed emission, which is 
the UV-heated gas. Two illustrative maps are shown for B335 
and L483mm in Fig.[^ 

Most sources show some excess emission on scales 

of 5"-10" or 1000-2000 AU at the average distance of 200 pc. 
The only exceptions are L1527 and OphIRS63. The emission 
is almost always aligned with one (12/24 sources) or both (4/24 
sources) outflow lobes. A few sources show widely distributed 
emission (6/24 sources). More Class 0 sources show excess 
emission along the direction of the outflow (11/13 sources) than 
Class Fs do (5/11) but this may bcsiS/N effect. 

The typical 6-5 line width is <1 km s“\ and so the 
emission is not part of the swept-up outflow gas as illustrated in 
more detail for the case of NGC 1333-IRAS 4A by jYildiz et al. 
( 2Q12| ). The only known mechanism to create this excess narrow 
emission is by UV photons generated from the protostellar ac¬ 
cretion process and subsequently escaping through the outflow 
cavities ( jSpaans et al.|1995| ). 

To estimate the effects of the UV radiation on these scales, it 
is first important to estimate t he teni perature of the gas compared 
with that of the dust. Figure A.14 shows model 3-2/6-5 
line ratios for a grid of kinetic temperatures and densities, with 
the observed values for each source overplotted at the 7.5" radius 
density. The inferred temperatures are in th e range of 30-80 K, 
consistent with the model predictions from j Visser et ak] ( |2012| ) 
on spatial scales of a few 1000 AU. For comparison, the typical 
dust envelope temperature at this distance is -15-25 K and thus 
the gas is heated to higher temperatures by more than a factor of 
2 . 

The mass of the UV-heated gas (component (ii)) is calculated 
on the basis of the residual after subtracting the model en¬ 
velope emission (component (/)) from the observed emis¬ 
sion. The mass is then calculated via the residual emission by 
assuming Fex^SO K and CO/H 2 =1.2x10“"^, where the value of 
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10 0 -10 10 0 -10 

Fig. 9. spectral maps in black overlaid with the model envelope 
spectra in red shown on the left panels. On the right panels, color maps 
of the UV heated gas distribution are shown. These are obtained by 
subtracting the model envelope emission from the observed spectra on a 
pixel-by-pixel basis. The sources are B335 (top) and L483mm (bottom). 
The axes show the offsets (Acr, A^) in arcsec. The color scale is in units 
of K km s"k 


Tpx is chosen because it is the median value for as reported 


t ex 

in 


Yildiz et ar] ( [2Q13| ) based on transitions from 2-1 up to 10- 


9. In order to compare UV-heated gas mass to the total outflow 
gas mass, the outflow mass is recalculated from 6-5 over 
the same (~30" x 30") area, using Tq^=75 K to be consistent 
with all other mass calculations. In Table the masses 
calculated for the envelope, UV-heated gas and outflow gas are 
tabulated. 

The mass of the UV-heated gas is typically a factor of 10 
to 100 times lower than the total envelope mass (Fig. [T^) and 
a factor of just a few up to 50 compared to the envelope mass 
within the 30"x30" region. There is no correlation with evolu¬ 
tion; i.e., the fraction of UV-heated gas compared to the total 
envelope mass does not change from Class 0 to Class I. Simi¬ 
larly there is no correlation between the mass of the UV-heated 
gas and the density at 1000 AU (Fig, p^), which may suggest 
that the emission is independent of density and thus the emission 
is thermalized. 

Compared with the outflow masses, the UV-heated gas 
masses (component (//)) are typically a few times higher, as also 
found for NGC 1333-IRAS 4A in |Yildiz et alT] ( |MT2| ). They fol¬ 
low a remarkably tight correlation with a Pearson correlation 
coefficient of 0.86 (3.1cr; Fig, p^). Furthermore, the fraction 
of UV-heated to envelope gas mass is constant as a function of 
bolometric luminosity at a median value of -0.03 (Fig.p^). The 
two outstanding high Muv/^env Class I sources are DK Cha and 
GSS30IRS1. 

Many protostellar envelopes show varying degrees of asym¬ 
metry and are not spheri cal; most striking is the flattened enve¬ 
lope surrounding LI 157 ( [Tobin et al.|2010| ). This asymmetry nat¬ 
urally introduces systematic uncertainties in the envelope model¬ 
ing which is then propagated through to the determination of the 
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Table 5. Comparison of envelope, photon-heated and outflow masses 
over the 30x30" area surrounding the central protostar. 


Source 

-^Envelope 

'^Envelope 

Muv" 

^outflow 


Total 

< 15" 

‘^CO 6-5 

‘2CO 6-5 

L1448mm 

9.0 

1.69 

4.0x10“^ 


NGC1333-1RAS 2A 

5.13 

0.67 

9.9x10-2 

8.1x10-2 

NGC1333-IRAS 4A 

5.59 

2.56 

8.2x10-2 

1.4x10-2 

NGC1333-IRAS 4B 

3.01 

2.60 

5.1x10-2 

2.2x10-2 

L1527 

0.92 

0.18 

2.3x10-2 

<5.5x10-'* 

CedllO-IRS4 

0.17 

0.07 

2.4x10-2 

3.3x10-'* 

IRAS 15398 

0.47 

0.28 

6.0x10-2 

5.5x10-'* 

L483mm 

4.4 

0.25 

3.9x10-2 

1.5x10-2 

Ser-SMMl 

16.13 

1.99 

3.0xl0-‘ 

7.4x10-2 

Ser-SMM4 

2.11 

3.17 

2.2xl0-‘ 

8.9x10-2 

Ser-SMM3 

3.21 

0.74 

2.2xl0-‘ 

6.1x10-2 

L723mm 

1.32 

0.67 

7.2x10-2 

1.9x10-2 

B335 

1.2 

0.79 

4.0x10-2 

9.8x10-'* 

L1489 

0.18 

0.04 

1.4x10-2 

1.5x10-'* 

TMRl 

0.23 

0.04 

1.6x10-2 

2.3x10-'* 

TMCIA 

0.22 

0.04 

7.6x10-2 

5.3x10-2 

TMCl 

0.2 

0.04 

7.7x10-2 

3.6x10-'* 

HH46-IRS 

4.36 

1.21 

2.5xl0-‘ 

1.2x10-2 

DK Cha 

0.82 

0.26 

2.2x10-2 

2.2x10-'* 

GSS30-1RS1 

0.6 

0.03 

3.1xl0-‘ 

2.5x10-2 

Elias 29 

0.3 

0.01 

7.8x10-2 

5.6x10-'* 

Oph-IRS63 

0.25 

0.12 

3.0x10-2 

<1.2x10-'* 

RN091 

0.45 

0.05 

1.1x10-2 

1.5x10-'* 

L1551-IRS5 

22.2 

0.23 

2.0x10-2 



Notes. All masses are given in M©. Total mass of the spherical en- 
velope inferred fr om the continuum radiative-transfer modeling ( |Kris-| 
jtensen et ^ 2012| ). As a, but limited to the mass within the 15" radius 
over which the UV-heated component is determined. UV photon- 
heated gas mass (component ( ii)) as calculated from the ^^CO 6-5 spec¬ 
tra over the mapped area after subtracting the modeled envelope emis¬ 
sion. Outflow mass calculated from the ^^CO 6-5 map for the same 
area as the ^^CO maps. All masses except the total envelope mass (a) 
are obtained over a ~30"x30" (Fig. |A.8|A.13j ) 


mass of the UV-heated gas. However, most envelopes are elon¬ 
gated perpendicular to the direction of the outflow (e.g., LI 157) 
whereas the residual ^^CO 6-5 emission is typically elongated 
along the outflow direction. Therefore we do not think that the 
use of spherical envelope models changes any of the conclusions 
regarding the effects of the UV-heated gas. 


4. Discussion 

4.1. Mechanical feedback 


Our results show that the outflow parameters inferred from the 
CO 6-5 data show the same trends with Lboi and evolutionary 
stage as found previously in the literature, but with stronger cor¬ 
relations than for the 3-2 data. Even though the same telescope 
and methods are used for all sources and the spatial resolution 
is high, there remains a scatter of at least an order of magnitude 
in the correlation between Fco and Lboi- This could point to the 
importance of “episodic accretion” as a resolution to the “lumi¬ 
nosity problem” ( Evans et al.|2009[[Dunham et al.|2010[|2013 ). 
Some Class 0 sources are very luminous, which is likely due to 
a current r apid burst in accretion which may happen every 10^- 
10"^ years ( Dunham et al.||2010 ). However, their location in the 
high state is not constant and would drop in the course of time, 
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Fig. 10. (a) UV heated gas mass is shown as a function of envelope mass (Menv), (b) density at 1000 AU (wh (1000 AU)), and (c) the outflow 
mass calculated from the 6-5 lines for the same region (Mobs(i 2 co 6 - 5 ))- (d) The fraction of the UV heated gas mass over envelope mass as a 
function of bolometric luminosity (Lboi). Figure p^d) has the same y-axis values as (a-c). 




Fig. 11. Correlation between Fco and the CO 14—13 and 18—17 fluxes 
obtained from Herschel-FACS. The integrated intensities are scaled to 
a common distance of 200 pc. The green lines shows the best-fit power- 
laws to the data and are simple least-squares fits. 


on timescales as fast as 10^ years ( Johnstone et al.||2013 ). The 
envelope mass, on the other hand, is independent of the current 
luminosity, and the stronger correlation of Fqo with Mgnv may 
simply reflect that more mass is swept up. 

Since the outflow force gives the integrated activity over 
the entire lifetime of a YSO, it is also interesting to compare 
this parameter with the currently shocked gas probed by the 
Herschel-FACS high-/ CO observations (/up>14). In Fig. | 
Fco is plotted against CO 14 -13 and CO 18-17 flu x es - 
and 940 K) obtained from [Karska et al.| ([2013]); jGoicoechea| 


et al. (|2012|); [Herczeg et aL| p012 ); Green et al.f ( |2013 ) and 


van Kempen et al.| (|2010al). there is a strong correlation with 
the CO 14-13 and CO 18-17 fluxes with Fco (r=~0.76 ~3.1o-; 
Fig. [TTJ. This correlation illustrates that although CO 18-17 
likely traces a different outflow component than CO 6-5, a com- 
ponent closer to the shock front ( [Santangelo et al.|2012t |Nisini| 
|et al.||2Ql'^ [Tafalla et al.||2Q13| ), the underlying driving mech¬ 
anism is the same. Furthermore, CO 18-17 emission is often 
extended along the outflow direction ( [Karska et al.||2013"] ) and 
clearly traces, spatially, a component related to that traced by 
CO 6-5. Although the excitation of CO 18-17 requires higher 
densities and temperatures (^crit~10^ cm“^; £’up~940 K) than 


CO 6-5 (^crit~10^ cm“^; £'up~120 K), CO 6-5 likely follows in 
the wake of the shocks traced by the higher-/ lines and therefore 
the excitation of both lines ultimately depend on the actual shock 
conditions. Testing this scenario requires velocity-resolved line 
profiles of high-/ lines such as CO 16-15 (Kristensen et al. in 
prep.). 


Fig. 12. Correlation between Fqo measured from the CO 6-5 (left) and 
CO 3-2 {right) data and the integrated intensity of the ground-state H 2 O 
lio-loi transition at 557 GH z. The integrated intens ities are scaled to a 
common distance of 200 pc ( [Kristensen et al.|2012j ). The correlation is 
strong for CO 6-5, r=0.90 with 4.1cr. The green lines shows the best-fit 
power-laws to the data and are simple least-squares fits. 


Another indication that the outflow force as measured from 
CO 6-5 is more closely linked to the currently shocked gas than 
that from 3-2 comes from comparing H 2 O and Fco- Water is 
one of the best shock trace rs, as shown most recently by sev- 
eral Herschel observations ( van Kempen et al.][2010b[ Lefloch 


et al. 20 101 [Kristensen et al.|2010[[Nisini et al.|2010t[Vasta et al. 


2012 Tafalla et al.|2013[ ). Kristensen et al. ( [2012 ) compared the 
integrated intensity of the H 2 O lio-loi transition at 557 GHz 
with the outflow forces presented in the literature. These ob¬ 
served line intensities are scaled by the square of the source 
distance to a common distance of Javg=200 pc. The literature 
values of the outflow force used in that paper were calculated 
using a variety of methods and data sets, and provided an inho¬ 
mogeneous sample. No correlation of H 2 O integrated intensity 
with Fco was found. Revisiting this comparison with the newly 
measured outflow forces in a consistent way reveals a correlation 
with the force measured from CO 3-2 data (r=0.78; 3.6cr) and 
a stronger correlation with the force measured from the CO 6-5 
data (r=0.90; 4Acr, Fig.[T^ (see also [Bjerkeli et al.|2012[ ). Thus, 
Fco as deduced from 6-5 can be used as a measure of the out¬ 
flow force of the shocked gas, rather than just the entrained, 
swept-up gas. 
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4.2. Radiative feedback 


The observational data demonstrate that 6-5 traces UV- 
heated gas, and that the UV-heated gas is predominantly found 
along the same direction as the outflow. The Ay is lower inside 
the outflow cavity, because the density is lower, and so UV radi¬ 
ation from the accretion can escape more easily along this direc¬ 


tion ( [Spaans et al.|1995[|Bruderer et al.|2009] Visser et al.|2012| ). 

If there are also external UV sources, the UV-heated gas could 
have a more isotropic component as well but this is not traced by 
our 6-5 data at the current S /N level except for the case 
of two sources in Ophiuchus, Elias 29 and GSS30IRS1. Narrow 
lines may be used instead at positions well away from the 


outflow cone, as illustrated by previous observations ( [van Kem- 
pen et al.|2009b| ). 


The estimated gas temperature of the UV-heated gas of 
30-50 K is likely a lower limit to the maximum temperature 
achiev ed b y this process. Mod el calculations by [Spaans et al. 


( 1995[ ) and [Visser et al.[p012[ ) show that the gas temperatures 


can reach values up to a few hundred K at 1000 AU radius in a 
narrow layer along the outflow cavity, depending on source char¬ 
acteristics. Gas temperatures >30 K are maintained out to 10"^ 
AU radius. Thus, in clustered environments such as NGC 1333 
or Ophiuchus, it is unlikely that the gas temperature ever drops 
down to 10 K because the protostars heat the gas radiatively. 
Gas and dust temperatures are clearly decoupled, with dust tem¬ 
peratures signiflcantly lower than the gas temperature, by about 
a factor of 2. Thus, estimates of the radiative feedback based 
on dust observations alone ( [Hatchell et al.|201^ [Sicilia-Aguilar[ 
et al.[[2013] ) likely underestimate both the temperature and the 
extent of the feedback. Indeed, the continuum emission, as ob¬ 
served with e.g. SCUBA at 450 and 850 yum, typically does not 
show extended structure along the direction of the outflows. 

The tight correlation between the mass of the UV-heated gas 
and the outflow mass, when measured over the same area, is puz¬ 
zling. Naively, one would expect the two properties to be unre¬ 
lated as they are caused by two different physical mechanisms, 
UV excitation and outflow entrainment. However, the cause of 
these two physical processes is linked, accretion and ejection 
(e.g. [Bontemps et al.|l 996} [Frank et al.|2014[ ). The UV photons 
are generated in the accretion shocks onto the protostar, and dur¬ 
ing this accretion process part of the material is ejected. Thus, 
higher-luminosity sources at a given envelope mass should show 
higher UV luminosities. It is not possible to verify this hypothe¬ 
sis directly as all of these sources are deeply embedded. A sec¬ 
ond component is required to efficiently UV-heat the surround¬ 
ing gas: an outflow cavity needs to be cleared out for the radia¬ 
tion to escape which requires the outflow to have been active for 
at least one dynamic time-scale. Thus, there may be good reason 
to expect a correlation between the masses of the UV-heated and 
outflowing gas, when measured over the same area. 


5. Conclusions 

In this paper, we present large-scale maps of 26 YSOs 
obtained with the APEX-CHAMP^ instrument (^^CO and 
6-5), together with the JCMT-HARP-B instrument (^^CO 
and 3-2). Our sample consists of deeply embedded Class 0 
sources as well as less deeply embedded Class I sources. With 
these high spatial and spectral resolution maps, we have studied 
the outflow activity of these two different evolutionary stages of 
YSOs in a consistent manner. All embedded sources show large- 
scale outflow activity that can be traced by the CO line wings, 
however their activity is reduced over the course of evolution 
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to the later evolutionary stages as indicated by the decrease of 
several outflow parameters, including the spatial extent of the 
outflow as seen in the ^^CO 6-5 maps. 

One of the crucial parameters, the outflow force, Fqo is 
quantified and correlations with other physical parameters are 
sought. In agreement with previous studies. Class 0 sources 
have higher outflow forces than Class I sources. Fqo is directly 
proportional to Mgnv and Moutfiow. showing that higher outflow 
forces are associated with higher envelope mass and outflow 
mass, as present in Class 0 sources. Comparing the outflow force 
as measured from CO 6-5 data to H 2 O observed with Herschel- 
HIFI and high-/ CO observed with Herschel-VACS reveals a 
correlation, suggesting that the outflow force from 6-5 is re¬ 
lated to current shock activity. This is in contrast with the outflow 
force measured from CO 3-2, where the correlation with water 
and the high-/ CO fluxes is weaker. 

The quiescent gas is traced by narrow (FWHM ~1 km s“^) 
^^CO 6-5 emission. For this purpose, maps are obtained in ^^CO 
6-5 transition for the sources ~1' region around the source posi¬ 
tion. Envelope emission is modeled via radiative transfer models 
and is subtracted from the observed ^^CO 6-5 emission. It is 
shown that an excess emission exists in most sources on scales 
of a 1000-2000 AU and this emission is caused by UV pho¬ 
tons generated from the protostellar accretion process and sub¬ 
sequently escaping through the outflow cavities. The fraction of 
the UV-heated gas compared to the total envelope mass does not 
change from Class 0 to Class I and there are no clear signs of 
evolutionary trends. 

UV heating is prominent along the outflow direction and this 
is a general observable trend. This directional preference sug¬ 
gests that the UV feedback on large scales is most important in 
the same regions as the outflows. The UV heating observed in 
6-5 is important on scales of <10"^ AU, i.e., not on cluster 
scales. Future models of core and disk fragmentation should take 
these effects into account. 
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Table A.l. Integration limits and contour levels. 





Blue Lobe" 



Red Lobe" 


CO 6-5 

CO 3- 

-2 

Source 

f^LSR 
[km s-i] 

Emax 
[km s“^] 

Yut,blue 

[km s-i] 

Yn,blue 

[km s“^] 

Emax 
[km s-i] 

t^in,red 

[km s“^] 

Y)ut,red 

[kms-1] 

Lowest Cntr* 
[Kkms-i] 

Step Size^ 

[K km s-i] 

Lowest Cntr^ 
[Kkms-i] 

Step Size^ 
[Kkms-i] 

L1448mm 

+5.2 

50.6 

-45.4 

2.0 

28.8 

6.6 

34.0 



10 

10 

NGC1333-IRAS 2A 

+7.7 

23.2 

-15.5 

7.0 

17.3 

10.5 

25.0 

25 

20 

15 

15 

NGC1333-IRAS 4A 

+7.0 

22.0 

-15.5 

7.0 

19.8 

9.2 

27.0 

20 

20 

20 

20 

NGC1333-IRAS 4B 

+7.1 

20.0 

-12.8 

4.0 

12.8 

9.2 

20.0 

20 

20 

20 

20 

L1527 

+5.9 

7.4 

-1.5 

4.5 

8.1 

7.0 

14.0 

7 

4 

5 

3 

CedllO-IRS4^ 

+4.2 

4.2 

0.0 

3.5 

3.8 

5.5 

8.0 

8 

5 

5 

3 

BHR7R 

-4.4 

15.5 

-20.0 

-6.0 

22.4 

-3.9 

18.0 

20 

20 



IRAS 15398 

+5.1 

8.6 

-3.5 

4.0 

9.9 

6.5 

15.0 

5 

5 

3 

3 

L483MM 

+5.2 

10.7 

-5.5 

3.5 

10.8 

6.0 

16.0 

8 

8 

5 

8 

Ser-SMMl 

+8.5 

19.0 

-10.5 

6.0 

22.5 

10.5 

31.0 

15 

20 

30 

25 

Ser-SMM4 

+8.0 

19.0 

-10.5 

6.0 

11.5 

10.5 

20.0 

15 

20 

30 

25 

Ser-SMM3 

+7.6 

22.0 

-13.5 

6.0 

13.5 

10.5 

22.0 

15 

20 

30 

25 

L723MM 

+11.2 

14.2 

-3.0 

9.0 

10.8 

12.0 

26.0 

15 

10 

5 

5 

B335 

+8.4 

8.9 

-0.5 

7.0 

8.6 

9.5 

17.0 

5 

5 



L1157 

+2.6 

12.4 

-9.8 

1.5 

17.4 

3.7 

20.0 



10 

20 

L1489 

+7.2 

13.7 

-6.5 

5.0 

7.8 

8.5 

15.0 

3 

2 

10 

10 

L1551-IRS5 

+6.2 

9.7 

-3.5 

4.5 

11.8 

7.5 

18.0 



10 

10 

TMRl 

+6.3 

7.8 

-1.5 

4.0 

3.7 

6.5 

10.0 

3 

3 

2 

2 

TMCIA 

+6.6 

17.1 

-10.5 

5.0 

5.4 

6.5 

10.0 

5 

5 

2 

5 

TMCl 

+5.2 

13.7 

-8.5 

4.0 

14.8 

6.2 

20.0 

4 

5 

3 

3 

HH46-IRS 

+5.2 

5.5 

-0.3 

10.0 

14.8 

12.2 

20.0 

20 

20 

10 

10 

DK Cha^ 

+3.1 

5.3 

-2.2 

1.5 

8.9 

4.3 

12.0 

5 

5 

5 

5 

GSS30-IRS1 

+3.5 

13.0 

-9.5 

1.5 

14.5 

5.5 

18.0 

20 

30 

15 

15 

Elias 29 

+4.3 

9.8 

-5.5 

1.5 

10.7 

7.0 

15.0 

15 

10 

7.5 

5.0 

Oph-IRS63 

+2.8 

11.3 

-8.5 

1.0 

3.2 

4.0 

6.0 

3 

1.5 

2 

1 

RN091 

+0.5 

15.0 

-14.5 

-1.0 

3.5 

1.0 

4.0 

3 

3 

3 

3 


Notes. Velocities are not corrected for inclination; i^LSR values are from |Yildiz et al. 
Contour levels are given in absolute intensities. Obtained from 6-5. 


( |2013| ). Velocity integration limits as shown in Fig. 
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Fig. A.l. CO 3-2 spectra with selected integration limits indicated, except for CedllOIRS4, BHR71, and DKCha where CO 6-5 was used. Each 
panel presents these limits for each source. The black spectrum at the bottom is taken from a clean position representative for the envelope emission. 
The blue spectrum at the middle is the representative spectrum from the blue outflow lobe, and red spectrum at the top is the representative spectrum 
from the red outflow lobe. Each panel shows flve vertical lines, these are Vlsr (black dashed line), Vout,biue (dot-dash blue line), Vin,biue (dashed 
blue line), Vin,red (dashed red line), and Vouued (dot-dash red line). 
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Fig. A.2. Each row contains contour and integrated intensity maps (in K km s of sources in 6-5 and 3-2. The contour levels and 
integration limits are given in Table |A.l] and integration limits shown in Fig. |A.l| The color images show all emission integrated from Vout,red to 
Vout,biue. including any minor cloud contribution. 
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Fig. A.3. Caption is same as for Fig. A.2 
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Fig. A.4. Caption is same as for Fig. A.2 
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Fig. A.5. 6-5 and 3-2 integrated intensity maps of the sources (in K km s ^). 
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Fig. A.6. Caption is same as for Fig. A.5 



Fig. A.7. Correlations between Fqo as measured from CO 3-2 and bolometric luminosity, envelope mass and outflow mass as determined from 
CO 3-2. 
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Fig. A.8. 7x7 pixel fully sampled maps are extracted toward the central positions of the sources in 6-5 (left) and in 6-5 (middle) 
transitions. The axes represent the equatorial offsets (Acr, A^) in arcsec. The main beam temperature intensity scale of each box are shown in 
the y-axes of the bottom-left box in Kelvins. The velocity range in each box is +8 km s"^ for the spectra, and +25 km s"^ for the 
spectra. The red lines in the left-hand panels are the 6-5 model line intensities for the passively heated envelope. The excess emission in the 
observations compared with these model profiles corresponds to the UV-heated gas and is shown as an image in the right p anel with the intensity 


scale in K km s ^. The middle and right panels contain the red and blue outflow lobes with the contour levels given in Table 
arrows in the right-hand panels show the direction of the outflow lobes. 
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Fig. A. 10. Same as Fig. A.8 
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Fig. A.ll. Same as Fig. A.8 
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Fig. A.12. Same as Fig. A.8 
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Fig. A.13. Same as Fig. 
blank. 
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6-5 transitions were not observed for L1448MM and L1551IRS5 in our observing campaign, therefore left 
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Fig. A.14. 3-2/6-5 intensity ratio as a function of density and 

gas temperature calculated via RADEX for A(^^CO)=1.5x10^"^ cm"^. 
Red markers indicate the observed intensity ratios for the central pixels 
for Class 0 sources whereas blue markers are for Class I sources. Both 
pixels are taken to be 15" diameter. The corresponding densities are the 
values at the 7.5" radi us found in the power-law envelope models of 
[Kristensen et al.| ( |2012| ). 












